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SECTION  1 
INTRODUCTION 


The  primary  goal  of  this  program  was  to  provide  a  simplified,  accurate 
method  of  neutron  dosimetry  for  radiobiological  studies,  so  that  quantitative 
assessments  of  radiobiological  effects  would  not  be  limited  by  the  accuracy  of 
the  dose-measurement  method  (the  cause),  but  rather  by  the  less  reproducible 
bio-effects  (the  effect).  The  probability  that  such  a  method  will  be  put  to 
use  is  a  direct  function  of  its  accuracy  and  an  inverse  function  of  its  costs 
in  turnaround  time,  manpower,  materials,  and  required  expertise  or 
sophistication.  The  method  as  presented  in  Section  4  below  incorporates  four 
neutron  activation  foils  with  natural  thresholds  extending  from  about  0.5  MeV 
to  8.7  MeV,  plus  cadmium-covered  and  bare  gold  foils  (0.414  x  10~6  MeV  and  0 
MeV  "thresholds").  It  utilizes  simple  hand-processing  of  the  data  and  is  the 
order  of  6%  accurate. 

Section  2  presents  the  "calibration"  of  the  ASTM  Standard  threshold-foil 
1-4 

spectrometry  method  in  the  NBS  Standard  Neutron  Field,  to  establish  its 
accuracy.  This  method  employs  the  SAND  II  unfolding  code,  and  was  used  in 
AFRRI  monkey  phantom  measurements  to  provide  a  basis  for  evaluating  the  simple 
5-foil  (plus  cadmium-covered  gold)  dosimetry  method  presented  in  Section  4. 

Section  3  presents  the  results  of  this  same  type  of  spectrum/dose 
measurement,  as  carried  out  at  several  points  inside  the  monkey  phantom.  The 
measurements  were  carried  out  both  with  the  standard  full  dosimetry  packet,  the 
"full  stack"  (i.e.,  the  baseline  measurement),  and  with  a  much  smaller, 
simplified  packet,  the  "short  stack."  The  accuracy  of  the  short  stack,  as 
compared  to  the  full  stack.  Is  demonstrated.  In  this  section  the  dosimetry 
results,  utilizing  the  SAND  II  spectrum  unfolding  code,  are  compared  to  the 
calculated  data  for  absolute  spectral  flux/dose. 

Section  4  presents  the  development  of  the  simplified  "short  stack"  method 
and  assesses  the  accuracy  of  this  straightforward,  hand-calculational  method  of 
threshold-foil  neutron  dosimetry.  The  prescription  for  carrying  out  accurate 
dosimetry  measurements  with  the  simplified  "short  stack"  method  is  also 
presented  here. 


SECTION  2 

CROSS-CALIBRATION  IN  NBS  STANDARD  NEUTRON  FIELD 


A  set  of  seven  (7)  neutron  activation  foils  was  irradiated  in  the  NBS 
Standard  Neutron  Field  in  order  to  cross-cal ibrate  the  threshold  foil  method  of 
neutron  spectrometry.  The  foils  were  counted  and  the  data  were  input  to  the 
SAND  II  unfolding  code  to  obtain  the  neutron  spectrum.  A  trial  spectrum  was 
also  input  to  SAND  II.  This  was  the  evaluated  U- 235  fission  neutron  spectrum 
discussed  below  (Figure  1). 

Table  1  presents  the  neutron  foil  activations,  as  well  as  corrections  in 
the  NBS  Standard  Neutron  Field  for  the  flux  activations  gradient  and  for 
self-scattering  by  each  foil. 

Figure  1  shows  the  SAND  II  output  spectrum  as  well  as  the  evaluated  U - 235 

fission  neutron  spectrum  (Designation:  XU5-5N1),  corresponding  to  a  fluence  of 
13  2 

7.786  x  10  n/cm  ,  +  1.7%.  The  corresponding  measured  value  utilizing  the 
ASTM  Standard  methodology1'4  yielded  a  fluence  of  7.98  x  1013  n/cm2,  which  is 
within  2.5%  of  the  NBS  value:  this  is  the  first  and  only  time  the  ASTM  Method 
used  here  has  been  so  evaluated. 

Thus,  the  accuracy  of  the  threshold-foil -spectrometry  method,  as  carried 
out  here,  appears  to  be  in  the  vicinity  of  2-3%  for  a  fission-type  spectrum. 
For  a  water-cooled,  water-moderated  reactor  spectrum,  a  1/E  tail  exists  below 
the  fission-like  spectrum.  However,  this  1/E  shape  is  pretty  much  a  canonical 
appendage;  i.e.,  it  is  known  to  rather  good  accuracy.  Therefore,  the  projected 
accuracy  for  a  water-moderated  reactor  is  probably  the  order  of  3-4%,  based  on 
the  NBS  calibration  results. 

It  should  be  pointed  out  that  while  the  SAND  II  "calibration"  produced  a 
value  of  neutron  fluence  within  2.5%  of  the  NBS  standard,  fluences  obtained 
from  some  of  the  individual  detectors  showed  variations  approaching  10%.  The 
strength  of  the  multiple-foil  approach  is  apparent  here  in  that  those  foils 
departing  furthest  from  the  average,  after  folding  in  each  one  with  the  trial 
spectrum,  are  weighted  the  least.  The  final  spectrum  was  obtained  from  the 
evaluated  spectrum  (used  as  input  to  SAND  II)  with  only  one  very  mild 
iteration.  The  spectral -shape  constraint,  as  well  as  the  de-emphasis  by  the 
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Figure  1.  ER1  Exposure  Room  Top  View,  Rhesus  Monkey  Phantom. 
ER2  is  same  but  without  Lead  Wall. 


code  of  those  foils  that  are  in  greatest  disagreement  in  terms  of  the 
least-squares  fit,  result  in  a  multi-foil  accuracy  considerably  greater  than 
that  implied  by  the  worst-foil  variation.  The  numerous  contributors  to  the 
errors  in  the  NBS  test  include  the  following. 

1.  Uncertainty  in  NBS  absolute  field  strength. 

2.  Uncertainty  in  the  flux-gradient  and  self-absorption  corrections. 

3.  Evaluated  fission-spectrum-shape  uncertainties. 

4.  Errors  Inherent  in  histogram-type  approximation  to  the  trial 
spectrum  (the  evaluated  spectrum). 

5.  Threshold-foil  cross  section  uncertainties. 

6.  Gamma-ray- 1 ine  measurement  uncertainties,  including  NBS-traceable 
calibration  gamma-ray  source  uncertainties,  fluorescent-yield 
uncertainties,  fission  branching  ratio  uncertainties  (for  fission 
foils),  activation  gamma-ray  line  decay  time  corrections,  and 
counting  statistics. 

Efforts  were  made  to  reduce  many  of  the  uncertainties  related  to  Item  6 
above.  A  large  improvement  was  achieved,  as  evidenced  by  the  fact  that  all  of 
the  i>AND  II  results,  those  for  the  NBS  foils  as  well  as  those  for  all  the  AFRRI 
phantom  spectra  presented  below,  required  only  one  or  two  mild  spectral -change 
perturbations;  a  large  improvement  over  earlier  measurements^  of  the  same  kind. 

In  summary,  an  accuracy  estimate  of  3-4%  for  the  AFRRI  monkey  phantom 
measurements  appears  to  be  quite  reasonable.  The  monkey  phantom  spectrum 
differs  from  the  NBS  U-235  fission  spectrum  in  that  the  former  has  a  1/E  tail 
plus  thermal -neutron  peak.  However,  the  1/E  tail  has  a  well-known  and 
accurately  calculable  shape,  and  only  contributes  about  11%  to  the  ER1, 
mid-thorax  dose.  The  thermal  peak  is  not  part  of  this  accuracy  evaluation,  as 
It  can  be  handled  both  experimentally  and  calculationally  as  an  isolated  case. 
(See  Section  6  for  details.) 


SECTION  3 

SAND  II  MEASUREMENTS  IN  PHANTOM 


This  section  presents  a  series  of  six  in-phantom  spectrum  measurements 
carried  out  with  the  short  stack  of  dosimetry  foils,  one  of  which  was  also 
carried  out  with  the  full  complement  of  foils.  This  will  be  referred  to  as  the 
full-stack  meeasurement  or  the  basepoint  measurement:  It  is  the  measurement 
which  serves  as  a  basepoint  for  the  simplified  5-foil  dosimetry  method 
presented  in  Section  4.  (See  Tables  2  and  3  for  the  full -stack  and  the  six 
short-stack  foil  data,  respectively.) 

The  accuracy  of  the  basepoint  measurement  was  established  in  Section  2 
above  where  all  but  the  (unnecessary)  low-energy  neutron-sensitive  foils  were 
used  for  measuring  a  hard  (unmoderated)  fission  spectrum.  In  this  section,  the 
accuracy  of  the  short  stack  is  evaluated  by  carrying  out  the  SAND  II  unfolding 
both  with  the  short  stack  and  the  full  stack  at  the  basepoint;  i.e.,  at  the 
mid-thorax  point  in  ER1  (Exposure  Room  1),  where  a  2"-thick  lead  gamma-ray 
shield  was  placed  between  the  reactor  and  the  phantom.  (Since  0"  to  6"  of  lead 
shielding  has  been  used  for  in-phantom  dosimetry  measurements,  the  ER1  exposure 
turned  out  to  be  a  good  choice  for  the  basepoint  measurement.  This  is 
especially  important  for  the  simplified,  hand-calculated  dosimetry  method 
developed  here  and  presented  in  Section  4  below  because  the  2"  lead  wall 
neutron  spectrum  is  intermediate  in  "softness"  between  the  0"  and  6" 
lead-shield  neutron  spectra.) 

Figure  2  shows  the  experimental  configuration  for  the  ER1  area,  where  a 
2"-thick  lead-wall  gamma-ray  shield  was  used.  The  ER2  configuration  was 
Identical,  except  that  no  lead  shield  was  used  there,  and  the  center  of  the 
phantom  was  100  cm  from  the  front  of  the  reactor  dimple  in  ER2,  instead  of  130 
cm. 

Figure  3  shows  a  cutaway  view,  through  the  rhesus  monkey  phantom  and 
perpendicular  to  the  reactor-phantom  center  line,  and  Figure  4  shows  the  simple 
phantom  and  chair  within  the  plywood  enclosure. 


Table  2.  Basepoint  Run  (All  Foils  Present)  Activations 

ER1:  Mid-Thorax. 


Foil* 

Enth(Mey) 

Ey(MeV) 

Mhr*1) 

NQX(sec 

157*u 

0+ 

0.412 

1.071-2++ 

0.95 

6.15-16 

«Mn 

0+ 

0.847 

2.687-1 

0.99 

3.66-16 

1.811 

2.687-1 

0.29 

235„ 

0+ 

0.537 

2.256-3 

0.26 

0.0626 

0.743 

4.126-2 

0.92 

0.0648 

8.93-11 

1.596 

1.723-2 

0.96 

0.0626 

237Np 

0.50 

0.743 

4.236-2 

0.92 

0.0568 

4.11-12 

1.596 

1.723-2 

0.96 

0.0548 

115Io 

1.00 

0.335 

1.540-1 

0.50 

— 

2.25-17 

238u 

1.45 

0.537 

2.256-3 

0.256 

0.0588 

0.743 

4.126-2 

0.92 

0.0509 

7.62-13 

1.596 

1.723-2 

0.96 

0.0588 

232Th 

1.75 

0.537 

2.256-3 

0.256 

0.0760 

1.78-13 

0.743 

4.126-2 

0.92 

0.0534 

54Fe 

2.20 

0.835 

9.532-5 

1.00 

— 

5.87-21 

58Ni 

2.90 

0.810 

4.051-4 

0.99 

— 

3.20-20 

2S 

6.30 

1.369 

4.621-2 

1.00 

— 

6.08-20 

27»1 

8.70 

1.369 

4.621-2 

1.00 

— 

3.02-20 

»°Zr 

14.00 

0.910 

6.641-3 

0.99 

... 

2.04-21 

*  All  foils  In  0.020"  cadmium 

**  N  X  In  sec’1  Is  given  for  the  non-fission  foils  only 
N°  is  given  for  the  fission  foils 

4  Since  the  foils  are  enclosed  in  cadmium,  the  effective  threshold  for 
these  foils  is  actually  4.25-7  MeV 

44  Read  as  1.071  x  10’* 


Table  3.  Hinl-Run  (6  Foils  Only)  Activations 
Used  for  Unfolding. 


1 

NoXtsec*1): 

ER1 

No 

X  (sec"*) :  ER2 

Foil* 

** 

Thorax 

Head 

Hip 

Thorax 

Head 

Hip 

6.15*16 

5.65-16 

6.98-16 

1.51-15 

1.32-15 

1.78-15 

1151„ 

2.25-17 

2.10-17 

2.86-17 

6.59-17 

5.57-17 

7.68-17 

54f. 

5.87-21 

5.71-21 

7.27-21 

1.97-20 

1.60-20 

2.08-20 

58n> 

3.20-20 

3.03-20 

3.81-20 

1.08-19 

8.81-20 

1.19-19 

MMg 

6.08-20 

5.82-20 

6.81-20 

2.21-19 

1.75-19 

2.43-19 

»*1 

3.02-20 

3.07-20 

3.52-20 

1.08-19 

8.99-20 

1.15-19 

All  foils  in  .020"  cadmium  cover 

These  activations  are  taken  from  the  standard  (all  foils  present)  run 


For  a  complete  list  of  activation  parameters,  cf,  either  Table  1  or  2 
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Figure  2. 


View  Perpendicular  to  Reactor-Phantom  Center  Line:  ER1  and  ER2 
Exposure  Room  Configurations. 


The  threshold-foil  activation  data  for  the  full  stack  of  dosimetry  foils, 
the  full  complement  of  (ASTM  Standard)  neutron  activation  foils,  is  presented 
in  Table  2  for  the  rhesus  monkey  phantom  mid-thorax  spectrometry/dosimetry 
exposure  in  ER1,  behind  the  2"  lead  shield;  i.e.,  the  basepoint  measurement. 
The  threshold-foil  data  for  the  six  (6)  short  stack  exposures  are  presented  in 
Table  3,  with  the  ER1  mid-thorax  data  being  common  to  the  Table  2  data.  All 
foils  presented  in  Tables  2  and  3  were  cadmium-covered,  with  a  watertight 
0.020"-thick  cadmium  enclosure.  The  bare  gold  and  manganese  foils  were  exposed 
In  a  separate  but  otherwise  identical  reactor  run,  and  the  activations  are  not 
presented  in  the  interest  of  brevity. 

The  thermal  group  is  group  37,  and  the  fluence/dose  determination  is 
achieved  with  a  simple  hand  calculation  utilizing  the  gold-  or  manganese-foil 
cadmium-difference  activations:  the  methodology  of  Sections  3  and  4  does  not 
apply  to  the  thermal -neutron  fluence/dose.  This  is  adequately  covered  in 
Section  6,  in  terms  of  the  goals  of  this  work. 

Table  4  presents  a  comparison  of  the  SAND  II  spectrum  unfolding  code 
output  for  the  basepoint  measurement  as  carried  out  with  the  full  stack 
(standard)  and  the  short  stack  (labeled  "mini").  The  calculated  spectrum/dose 
data  are  also  presented  in  Tables  4,  5  and  6  for  calculations  carried  out  in 
the  basepoint  location  (ER1,  mid-thorax),  at  center  head  in  ER1,  and  for  both 
mid-thorax  and  center-head  locations  in  ER2.  The  comparison  of  measured  and 
calculated  data  is  presented  in  Section  3.2  below  and  includes  the  experimental 
normalization  of  the  calculational  data,  a  vital  point  when  considering  that 
one  would  otherwise  have  to  rely  on  chamber  dosimeters,  whose  accuracy  is,  in 
effect,  being  evaluated  by  these  more  sophisticated  spectrometry/dosimetry 
methods  presented  here. 

The  data  of  Table  4  show  an  agreement  of  2%  in  total  fluence  between  the 
full -stack  and  the  short-stack  SAND  II  code  outputs.  The  doses  agree  to  within 
about  IX.  Application  of  a  quicker,  analytical  technique  utilizing  these  same 
threshold-foil  activation  data  and  the  accuracy  of  that  method  are  presented  in 
Section  4  below.  That  accuracy  assessment  is,  of  course,  based  directly  on  the 
use  of  the  SAND  II  results  (Tables  4,  5  and  6),  and  on  the  adequacy  of  these 
baseline  data. 


I  4-ftft  y.i'l-l-t  -ft  -  ft  1  ■  ft  ft'a 


Table  4.  Comparison  of  Basepoint  Run  {All  Foils  Present)  with  Mini-Run  (6  Foils  Only) 


for  ER1,  Mid-Thorax  •  6roup  Fluences  (n/cnr 


Min),  Doses  (Rad  Tissue) 


Upper 


Calculated 


Basepoint 


YY 

Group  # 

Energy 

(MeV) 

*Calc 

DCalc 

^Meas 

°Meas 

*Meas 

^Meas 

1 

19.6 

6.64+3 

4.89-5 

7.42+3 

5.46-5 

3.70+3 

2.72-5 

2 

16.9 

3.60+4 

2.49-4 

5.21+4 

3.62-4 

2.64+4 

1.83-4 

3 

14.9 

3.10+4 

2.07-4 

4.19+4 

2.80-4 

2.15+4 

1.44-4 

4 

14.2 

1.98+4 

1.30-4 

2.95+4 

1.94-4 

1.53+4 

1.00-4 

v 

5 

13.8 

1.42+5 

9.02-4 

1.73+5 

1.10-3 

9.35+4 

5.95-4 

5 

6 

12.8 

1.30+5 

7.95-4 

1.55+5 

9.49-4 

9.66+4 

5.91-4 

» 

7 

12.2 

4.99+5 

3.03-3 

3.52+5 

2.14-3 

3.47+5 

2.11-3 

8 

11.1 

8.55+5 

4.90-3 

6.45+5 

3.70-3 

6.43+5 

3.69-3 

.  * 

1  'll 

9 

10.0 

1.51+6 

8.35-3 

1.10+6 

6.06-3 

1.09+6 

6.01-3 

i 

10 

9.05 

2.37+6 

1.25-2 

1.79+6 

9.43-3 

1 .78*6 

9.37-3 

11 

8.19 

4.37+6 

2.24-2 

3.23+6 

1.65-2 

3.20+6 

1.64-2 

•  V 

12 

7.41 

1.23+7 

6.08-2 

9.05+6 

4.48-2 

8.99+6 

4.45-2 

13 

6.38 

3.32+7 

1.53-1 

2.69+7 

1.24-1 

2.74+7 

1.26-1 

-•'4 

14 

4.97 

8.02+6 

3.51-2 

6.89+6 

3.02-2 

7.07+6 

3.09-2 

*v 

15 

4.72 

2.86+7 

1.23-1 

2.35+7 

1.01-1 

2.41+7 

1.04-1 

- '  a . 

16 

4.07 

6.75+7 

2.75-1 

5.87+7 

2.39-1 

6.07+7 

2.47-1 

$ 

17 

3.01 

7.74+7 

2.66-1 

6.39+7 

2.20-1 

6.64+7 

2.28-1 

18 

2.39 

1.34+7 

4.22-2 

1.10+7 

3.47-2 

1.14+7 

3.60-2 

19 

2.31 

8.19+7 

2.54-1 

7.03+7 

2.19-1 

7.38+7 

2.29-1 

■  *1 

20 

1.83 

2.02+8 

5.41-1 

1.61+8 

4.31-1 

1.66+8 

4.46-1 

t  / 

21 

1.11 

2.30+8 

4.70-1 

1.90+8 

3.90-1 

1.90+8 

3.89-1 

rV» 

22 

.550 

1.97+8 

2.51-1 

1.71+8 

2.17-1 

1.61+8 

2.04-1 

ft 

23 

.158 

4.57+7 

3.57-2 

3.60+7 

2.81-2 

3.13+7 

2.44-2 

24 

.111 

9.51+7 

5.17-2 

7.16+7 

3.90-2 

6.31+7 

3.43-2 

P 

25 

.0525 

8.11+7 

2.49-2 

6.04+7 

1.85-2 

5.32+7 

1.63-2 

*’  F 

26 

.0248 

1.05+7 

2.18-3 

8.23+6 

1.72-3 

7.22+6 

1.51-3 

27 

.0219 

7.93+7 

1.14-2 

5.60+7 

8.01-3 

4.96+7 

7.11-3 

f,'! 

28 

.0103 

9.04+7 

5.52-3 

6.93+7 

4.23-3 

6.22+7 

3.80-3 

29 

3.35-3 

8.49+7 

1.83-3 

6.49+7 

1.40-3 

5.85+7 

1.26-3 

‘1**> 

30 

1.23-3 

8.50+7 

7.66-4 

6.04+7 

5.45-4 

5.41+7 

4.88-4 

ri 

31 

5.83-4 

1.72+8 

5.36-4 

1.24+8 

3.87-4 

1.10+8 

3.43-4 

1, 

32 

1.01-4 

1.20+8 

1.46-4 

8.52+7 

1.03-4 

7.75+7 

9.40-5 

33 

2.40-5 

9.66+7 

1.22-4 

7.24+7 

9.12-5 

6.49+7 

8.20-5 

ft 

34 

1.07-6 

1.51+8 

2.96-4 

9.91+7 

1.94-4 

9.12+7 

1.78-4 

35 

3.06-6 

8.83+7 

2.93-4 

6.88+7 

2.28-4 

6.08+7 

2.02-4 

* 

36 

1.13-6 

9.59+7 

5.22-4 

5.76+7 

3.13-4 

5.10+7 

2.78-4 

>  i 

37 

4.14-7 

4.34+9 

9.72-2 

2.84+9 

6.37-2 

2.84+9 

6.37-2 

ji 

/«■ 

V* 

Total 

6.60+9 

2.76+0 

4.57+9 

2.25+0 

4.48+9 

2.28+0 

Vi'V.v.'Si 
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3.1  SPECTRAL  SHAPES  . 


The  SAND  II  unfolding  results  of  the  baseline  measurements  at  the 
basepoint  (ER1,  mid-thorax)  are  shown  for  the  standard  dosimetry  foil  pack  in 
Figure  5  and  for  the  short  stack  in  Figure  6.  They  are  compared  to  the  direct 
calculational  outputs:  these  calculational  data  were  normalized  to  some 
tissue-equivalent  (TE)  dosimetry  measurements  taken  at  some  other  (distant) 
time,  and  the  agreement/disagreement  with  the  absolute  measurement  is  therefore 
misleading.  A  valid  renormalization  of  the  calculated  data  was  done  utilizing 
nickel  foils  placed  on  the  front  of  the  reactor  dimple  (see  Figure  2),  and  a 
discussion  of  calculational  predictions  in  magnitude  is  deferred  to  Section 
3.2.  The  full-stack  and  short-stack  data  are  in  excellent  agreement  with  one 
another  in  shape  as  well  as  magnitude.  The  agreement  in  shape  with  the 
calculated  (trial)  spectrum  is  also  excellent,  as  evidenced  by  the  fact  that 
the  SAND  II  code  carried  out  its  very  minimum  of  only  one  iteration  to  obtain 
the  final  (measured)  spectrum  that  is  in  accord  with  the  threshold-foil  data. 

Figure  7  compares  the  absolute  measured  data  with  the  calculated  data,  as 
normalized  to  some  TE  dosimeter  measurement  carried  out  at  a  different  time. 
The  agreement  in  spectral  shape  is  indeed  very  good,  as  it  is  for  the 
measurement  in  the  rhesus  monkey  phantom  hip  in  ER1.  These  data  for  the  hip 
location,  shown  in  Figure  8,  utilized  the  neutron  spectrum  calculated  at 
mid-thorax  in  ER1.  The  calculated  spectrum  was  arbitrarily  fitted  to  the 
absolute  spectral  flux  measurement,  in  Figure  8,  to  show  a  good  shape 
comparison. 

Figures  9,  10  and  11  show  in-phantom  spectral  data  similar  to  that  of 
Figures  6,  7  and  8,  respectively,  but  for  the  free-field  case  (i.e.,  no  lead 
shield)  in  ER2.  Again,  the  measured  spectral  shapes  are  in  good  agreement  with 
the  calculated  shapes. 

In  Figure  12,  a  comparison  in  spectral  shape  is  shown  between  ER1  and  ER2 
spectra,  mid-thorax.  Note  that  the  spectral  shapes  are  practically  identical 
below  about  0.4  MeV.  This  is  the  "canonical  slowing-down  spectrum."  As 
normalized  in  the  slowing-down  region,  the  ER1  spectrum  above  0.4  MeV  is  seen 


UIUI  •  M>l  ->sd  A9W-  WO/U  ASW  ‘(3 )<P*  3 
6 


22 


.VA 


to  be  appreciably  softer  due  to  the  lead  wall.  Inelastic  scattering  in  the 
lead  shield  depletes  this  high-energy  region  from  0.4  MeV  on  up. 

3.2  COMPARISONS  OF  AFRRI  ZOOLOGICAL  PHANTOM  MEASUREMENTS 

WITH  CALCULATIONS  NORMALIZED  TO  NICKEL  FOIL  MONITORS  . 

Neutron  fluences  have  been  measured  and  calculated  for  mid-head  and 
mid-thorax  locations  in  a  simple  zoological  phantom  exposed  at  the  AFRRI  TR1GA 
reactor  facility.  The  calculated  data  have  been  produced  by  Kaul  et  al,6  per 
unit  free  field  KERMA  and  were  first  normalized  to  free  field  KERMA  measured  at 
the  intended  experiment  location,  100  cm  from  the  reactor  centerline  as 
follows: 

Exposure  Room  1  with  2"  lead  wall 

7.8  Rad(tis)  per  Kw-Min 

Exposure  Room  2  free  field 

10.8  Rad(tis)  per  Kw-Min 

The  experimental  measurements  in  ER1  were  performed  at  130  cm  from  the  core 
centerline  rather  than  100  cm.  Thus,  the  reported  calculations  had 
consequently  been  adjusted  by  Verbinski-Cassapakis  by  dividing  by  1.617,  a 
factor  determined  from  measurements  made  at  similar  reactors.  The 

calculational  results  were  substantially  greater  than  the  measurements, 
differing  by  factors  ranging  from  20  to  100%,  which  brought  to  question  the 
practice  of  normalizing  the  calculations  to  standard  neutron  dosimeters  and 
comparing  the  results  to  the  much  more  sophisticated  dosimetry  measurements 
presented  here. 

Because  these  calculation-measurement  discrepancies  are  much  larger  than 
those  reported  in  a  similar  comparison  reported  by  Kaul  et  al,  it  was  decided 
to  reexamine  the  measurement-calculation  comparison  based  on  a  normalization  to 
nickel  foil  activation  (Ni -58  (n,p)  Co-58)  data  taken  at  the  same  time  as  the 
phantom  measurements.  The  nickel  foils  were  taped  to  the  aluminum  "dimple" 
which  protrudes  Into  each  exposure  room  and  into  which  the  reactor  is  moved  at 
the  time  of  each  experiment  (see  Figure  2).  They  were  positioned  at  the 
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reactor  midplane  along  a  direct  line  to  the  phantom.  The  nickel  foil  results, 
in  terms  of  disintegrations  per  second  per  Kw-Min  per  Ni-58  atom,  at  exposure 
time,  are  as  follows: 


ERKW/Lead) 

5.33*10‘22 


ER2(M/Q  Lead) 

,-22 


6.11*10 


These  values  are  converted  to  atoms  of  Co-58  per  Kw-Min  per  atom  of  Ni-58  based 


on  a  Co-58  disintegration  rate  of  1.125*10 
,-15 


-7 


4.74*10 


5.43*10 


per  sec. ,  as  follows: 
15 


The  free  field  fluences  on  which  the  phantom  calculations  are  based  have  been 

7 


performed  using  ANISN,  a  one-dimensional  discrete  ordinates  code,  in  a 
spherical  geometry  model  which  includes  the  reactor  and  all  intervening  detail. 
The  calculated  neutron  scalar  fluences  at  the  nickel  foil  locations  have  been 
used  together  with  Vitamin-E/ENDF/B-V  activation  cross  sections8  to  calculate 
the  numbers  of  Co-58  nuclei  per  Ni-58  nucleus  per  source  neutron  which  are  as 
follows: 


1.76*10 


•30 


1.50*10 


-30 


Thus,  the  numbers  of  source  neutrons  per  Kw-min  required  to  calculationally 
reproduce  the  measured  activation  in  each  room  are: 


2.69*10 


15 


3.62*10 


15 


Using  these  values,  the  neutron  KERMA  In  units  of  Rad(tis)  per  Kw-Min  at  100  cm 
from  the  reactor  calculated  using  one-dimensional  methods  are: 


6.13 


9.27 


which  are,  respectively,  78.6%  and  85.8%  of  the  AFRRI  monitor  values  taken  at 
the  same  locations  but  on  a  different  day. 
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The  nickel  foil-based  source  normalizations  were  used  to  obtain  new 
phantom  fluences  for  comparison  with  the  experimental  values.  The  renormalized 
calculational  results  and  the  experimental  results  are  presented  in  Table  7. 
The  agreement  at  the  mid-thorax  location  for  neutron  energies  greater  than  10 
KeV  Is  10%  or  better.  Because  the  mid-head  location  is  actually  out  of  the 
reactor  mid-plane  It  Is  expected  that  the  calculated  values  will  show  an 
Increase  relative  to  those  of  the  measurements  (unless  the  calculations  are 
corrected  for  the  larger  distance  of  mid-head  to  reactor  center).  This  is  the 
case  for  Exposure  Room  2  in  the  bare  configuration.  However,  it  is  not  the 
case  in  Room  1,  which  has  the  lead  shield.  The  lead  appears  to  have  the  effect 
of  creating  a  more  uniform  exposure  over  the  phantom.  For  energies  below  10 
KeV  the  agreement  between  calculation  and  experiment  is  no  better  than  10%  and 
is  generally  much  worse.  This  is  most  likely  the  result  of  using  an  overly 
simplified  exposure  field  model. 

3.2.1  Summary  . 

Good  agreement  has  been  obtained  between  calculated  and  measured  fast 
fluence  and  KERMA  values  inside  a  zoological  phantom  exposed  at  the  AFRR1  TR1GA 
reactor.  This  agreement  has  been  obtained  on  the  basis  of  normalization  to 
nickel  activation  foil  monitor  data  taken  on  the  face  of  the  reactor 
containment.  Earlier  comparisons  in  which  calculations  were  normalized  to 
paired  chamber  monitor  data  taken  at  the  nominal  experiment  location  yielded 
much  poorer  agreement.  It  is  expected  that  better  agreement  can  be  obtained  at 
low  neutron  energies  and  at  locations  out  of  the  reactor  midplane  by  using  free 
field  fluences  obtained  by  two-dimensional  methods,  as  contrasted  to  the 
one-dimensional  calculations  presented  here. 


Tabl*  7 .  Four-Group  Co-par Iso*  of  Calculated  Neutron  Fluence  and  Dose  with  Baseline  Measurements 
Values  per  Kw-«1n  at  130  c-  (ER1)  and  100  oa  (ER?)  fro-  Reactor  Centerline  . 


L*  fc*l*  (tlMl'all'rf  VUi'tli*  Jft*  JjVjVt*. 


-“.V,  /.  V.  r.  v>,  r.  r,  <,  w,  r.y. 
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SECTION  4 

THE  5-FOIL  NEUTRON  DOSIMETRY  METHOD 


This  section  presents  a  simple  hand-calculational  (perturbation)  method  of 
threshold-foil  dosimetry  in  which  5  threshold-activation  foils  are  used.  (The 
SAND  II  code  Is  not  required;  neither  is  there  need  for  a  trial  spectrum. 
However,  the  "Incremental  doses"  assigned  to  each  foil  (Table  8)  were  derived 
from  the  SAND  II  basepoint  measurement.)  One  foil,  the  nickel  foil,  is  used  to 
derive  the  first-order  dosimetry  value,  D°(Ni).  This  Is  simply  equal  to  D  , 
the  dose  of  the  baseline  run  (i.e.,  the  "standard"  run  in  Table  4),  multiplied 
by  the  ratio  of  the  nickel -foil  activations. 

D° (Ni )  -  DB  R(Ni)  4.1 

where  R(Ni)  -  (NQX)X  (Ni  )/(NQX)B  (Ni ) . 

Here,  (NqA)x  refers  to  the  dosimetry  measurement  under  consideration  and 
(NqX)B  refers  to  the  baseline-run  value  (the  value  for  mid-thorax,  ER1,  in 
Table  3). 

Now,  If  there  were  no  change  of  spectral  shape  between  the  baseline  run 
and  run  x,  Equation  4.1  would  give  an  accurate  measure  of  the  dose:  The  nickel 
foil  has  proven  to  be,  by  far,  the  most  reliable  foil  in  terms  of  reproducible 
results. 

The  nickel  foil  by  itself  cannot,  of  course,  correct  for  a  change  of 
spectral  shape.  However,  the  foils  with  thresholds  at  energies  different  from 
the  nickel  threshold  of  about  3  MeV  do  provide  change-of-spectral -shape  data. 
This  is  done  by  means  of  the  following  "shape-perturbation"  correction. 


Here,  the  A.  refers  to  the  activation  ratios  of  the  i-th  foil  (i.e.,  Au,  In,  Fe 
and  Mg)  for  the  unknown  and  the  baseline  run,  divided  by  R ( Ni ) :  i.e., 

A(Au)  -  R(Au)/R(Ni ) 

where 

R(Au)  -  (N0X)X(Au)/(N0X)B(Au). 

The  various  refer  to  the  fraction  of  the  dose  "sensed"  by  the  i-th 
foil,  in  the  case  of  the  baseline  run,  excluding  the  nickel  foil.  The  F^ 
values  were  determined  from  the  dose  values  D^-j-q(R)  given  in  Table  8  for  the 
"standard"  or  baseline  run;  i.e.,  the  ER1,  mid-thorax,  full-stack  spectrum. 
The  corresponding  energy  range  (region  of  maximum  influence)  for  each  foil  is 
also  presented  in  Table  8.  Table  9  presents  the  value  of  F.  for  each  of  the  4 
"shaping-sensing"  foils. 

4.1  APPLICATION  OF  THE  SIMPLE  5-FOIL  DOSIMETRY  METHOD. 

The  foil -activation  values  presented  in  Table  3,  for  all  six  (6)  spectrum 
measurements,  were  utilized  to  provide  the  doses  in  phantom  by  means  of 
Equations  4.1  and  4.2  above.  The  dose-fractions  F ^ ,  Equation  4.2,  are  given  in 
Table  9  for  each  "shape-correction"  foil  i  (i.e.,  Au,  In,  Fe  and  Mg). 

The  values  of  A.  (i.e.,  R(Mg)/R(Ni))  are  presented  in  Table  9,  along  with 
D°(Ni),  the  first  order  nickel-foil  dose  determination,  and  Dx,  the  second 
order  determination  in  which  a  change-of-spectral -shape  correction  has  been 
applied.  Note  how  the  A^  systematically  Increase  with  neutron  energy  for  the 
ER2  spectra:  As  shown  In  Figure  11,  these  spectra  are  appreciably  harder  than 
the  ER1  spectra,  due  to  the  down-scattering  (in  energy)  of  neutrons  by  the 
2-inch  thick  lead  wall  in  ER1. 

The  values  of  Dx  are  presented  in  Table  9,  along  with  Ds,  the  values 
derived  from  the  SAND  II  unfolding-code  spectra.  The  last  row  presents  the 
percent  variation  of  Dx  from  Ds.  The  standard  deviation  for  the  ER1  (excluding 
mid-thorax)  and  ER2  data  are  6%. 
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Table  8. 

Standard 

* 

Run  Parameters 

for  Short-Stack  Dosimetry  at  AFRRI 

• 

Foil 

Represented 

Energy  Range 

EL(MeV)  Eu(MeV) 

(No^STD^sec  ^ 

N+ 

n  STD 

dstd(K) 

197Au 

4.14-7 

0.55 

6.68-16 

2.09+4 

1460.2 

115ln 

0.55 

2.31 

2.25-17 

7.03+2 

4551.4 

54Fe 

2.31 

3.01 

5.87-21 

1.83-1 

1114.6 

58Ni 

3.01 

6.38 

3.20-20 

1 

2142.9 

24Mg 

6.38 

19.6 

6.08-20 

1.90+0 

375.4 

* 

The  phantom  mid-thorax  exposure  behind  2"  lead  at  ER1  is  taken  to  be 
the  standard  run. 

N+stq  is  the  Ni-normalized  foil  activation. 

++  ( R)  are  doses  for  the  4187.5  kW  •  Min  exposure,  between  E^  and  Ey. 

Columns  2  and  3,  taken  from  the  basepoint  ("standard")  measurement,  Table  4, 
but  with  ICRP  rather  than  ICRU  dose  conversions.  The  ICRP  vs.  ICRU  differences 
are  irrelevant  here. 


Table  9.  Dos*  Oeteratnetlons  via  5-Foll  Ooslaetry  Stack  for  CpIUeraal  Neutrons,  per  kN  -  HU, 
Coapared  to  from  SANO  II  Unfolding  «lth  S  Foils.  Doses  In  Rad  (Tlssee)  per  kM>  Nln 
TRo  D#  and  0*  Values  are  Hand-Calculated  Values  (See  Equations  4.1  and  4.2). 


The  5-foil  stacks  for  the  mid-thorax  spectrum,  ER1,  were  taken  from  the 
baseline-run  data  presented  in  Table  2.  When  the  SAND  II  code  was  run  with  the 
5-foil  stack,  the  resulting  epi -thermal  neutron  dose  was  about  1%  higher.  This 
is  somewhat  representative  of  the  error  resulting  from  utilizing  a  small 
fraction  of  the  full -foil -stack  data.  This  is  to  be  compared  with  the  6% 
standard  deviation  resulting  from  utilizing  the  same  five  foils,  but  without 
SAND  II. 


SECTION  5 

COMPARISON  OF  SANO  II  (5-FOIL  STACK)  AND  THE 
SIMPLE  5-FOIL  DOSIMETRY  METHOD  PRESENTED  IN  SECTION  4 

The  Improvement,  due  to  the  use  of  SAND  II,  derives  from  making  full  use 
of  the  a  priori  knowledge  of  the  spectral  shape.  This  a  priori  knowledge  is 
Input  to  SAND  II  In  the  form  of  the  trial  spectrum;  l.e.,  the  spectrum  from  the 
best  available  calculation.  Thus,  SAND  II  also  represents  a  perturbation 
method  of  neutron  spectrometry/dosimetry.  The  best-known  (input)  spectrum  is 
perturbed  In  shape  and  absolutely  adjusted  In  magnitude  to  agree  with  the 
foil -activation  data. 

The  5-foil  dosimetry  method  presented  in  Section  4  can  be  thought  of  as  a 
perturbation  method  as  well  in  that  it  corrects  for  variations  in  spectral 
shape  via  an  approximate  5-group  ("histogram")  approach.  By  contrast,  SAND  II 
makes  full  utilization  of  the  a  priori  knowledge  of  the  spectral  shape. 
However,  the  estimated  6%  accuracy  for  the  5-foil  5-group  "histogram"  method 
may  be  adequately  accurate,  compared  to  the  3-4%  estimated  accuracy  for  the 
maxi-foil  dosimetry  stack  used  with  SAND  II.  Consider  further  that  no  accurate 
a  priori  knowledge  of  spectral  shape  is  needed  for.  the  5-foil  5-group 
"histogram"  method  of  Section  4,  making  it  much  easier  to  apply. 

The  estimated  6%  accuracy  is  to  be  contrasted  to  the  inaccuracies  implied 
in  Section  3.2,  where  the  calculational  results®  were  first  obtained  via 
normalization  to  measured  free  field  KERMA  (measured  with  TE  dosimeters). 
Discrepancies  ranging  from  20  to  100%  were  reduced  to  the  order  of  10%, 
according  to  Section  3.2.  This  section,  as  well  as  Section  3.2.1,  was  taken 
directly  from  a  communication  by  Dean  Kaul  (Letter,  August  22,  1985). 


SECTION  6 

THERMAL  (SUB-CADMIUM)  NEUTRON  DOSE 

The  thermal -neutron  dose  is  derived  from  gold  foil  data,  utilizing  the 
cadmium-difference  method.  The  gold-foil  cross  section  is  taken  as  78.54 
barns,  which  is  the  cross  section  derived  from  a  Maxwellian  thermal -neutron 
peak  plus  1/E  tail  extending  out  to  0.414  eV.  This  cross  section  value  is  more 
accurate  than  the  value  of  98.8  barns,  for  2200  m/sec  neutrons. 

A  value  of  2.24  x  10***  rads  per  unit  neutron  fluence  is  the  proper  dose 
conversion  factor  (ICRU  muscle  KERMA). 

It  must  be  pointed  out  that  the  bare  and  cadmium-covered  gold  (or 
manganese  or  cobalt)  foils  cannot  be  exposed  one  behind  the  other,  or  even 
side-by-side,  because  of  the  cadmium-produced  local  flux  depression:  If 
side-by-side,  but  well  separated  (say  by  3  or  4  cadmium  "diameters"),  then  it 
is  incumbent  on  the  experimenter  to  prove  that  both  are  in  identical  epicadmium 
neutron-flux  fields.  In  the  measurements  presented  here,  the  bare  foils  were 
exposed  in  separate  reactor  runs,  with  run-to-run  normalization  carried  out 
with  nickel  activation  foils  taped  to  the  reactor  entrance  window  (on  the  front 
of  the  "dimple"  in  Figure  2). 

These  nickel  monitor  foils  proved  to  be  a  fortunate  choice  because  they 
provided  a  basepoint  for  the  "reactor  power"  or  "reactor  source  term" 
renormalization.  This  is  an  otherwise  formidable  task,  requiring  flux  plots 
throughout  the  reactor,  an  exact  knowledge  of  day-to-day  secondary  standard 
(fission  chamber,  say)  calibration,  of  the  thickness  of  water  between  reactor 
and  exit  window  to  ER1/ER2,  etc.  (Prior  to  this,  the  calculation  was 
renormalized  to  some  dosimeters  used  at  a  much  different  time.  This  proved 
unfortunate,  as  can  be  seen  from  the  results  of  Table  4,  for  example. 
Renormalization  to  the  nickel  foils  on  the  dimple  provided  a  meaningful 
comparison  of  measurement  and  calculation,  which  would  have  otherwise  proven 
untenable. ) 
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